The disordered-free-moment phase: a low-field disordered state in spin-gap 

antiferromagnets with site dilution 
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Site dilution of spin-gapped antiferromagnets leads to localized free moments, which can order 
antiferromagnetically in two and higher dimensions. Here we show how a weak magnetic field drives 
this order-by-disorder state into a novel disordered-free-moment phase, characterized by the forma- 
tion of local singlets between neighboring moments and by localized moments aligned antiparallel 
to the field. This disordered phase is characterized by the absence of a gap, as it is the case in a 
Bose glass. The associated field-driven quantum phase transition is consistent with the universality 
of a superfluid-to-Bose-glass transition. The robustness of the disordered-free-moment phase and 
its prominent features, in particular a series of psetido-plateaus in the magnetization curve, makes 
it accessible and relevant to experiments. 
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Valence bond solids (VBS) in spin-gapped antiferro- 
magnets represent some of the most fundamental exam- 
ples of quantum-disordered states in condensed matter 
systems. The nature of such states is by now well under- 
stood theoretically and has been extensively verified ex- 
perimentally. A variety of mechanisms, such as increased 
strength of certain bonds in the magnetic Hamiltonian 
[3, [|| , magnetic frustration 0] , and the Haldane mecha- 
nism 0] , can render classical Neel order unstable towards 
the formation of local singlets, which arrange themselves 
into a VBS. Evidence for such phases has been found 
in a large number of magnetic compounds, ranging from 
Haldane chains to spin ladders (|, to weakly coupled 
dimer systems |7|, |8( . A central focus of theoretical and 
experimental investigations has been the effect of doping 
on such states. In particular, it was soon realized the- 
oretically, and observed experimentally [ll], H3 |. 
that doping a VBS with static, non-magnetic impurities 
leads to the intriguing phenomenon of order-by-disorder 
(OBD): free S — 1/2 magnetic moments appear close 
to the impurity sites and interact effectively via a long- 
range network of unfrustrated (albeit random) couplings, 
which decay exponentially with the inter-moment dis- 
tance. These interactions, although weak, are sufficient 
for the free moments (FMs) to order antiferromagneti- 
cally at experimentally relevant temperatures 

Given the random nature of the inter-moment cou- 
plings, the OBD state induced by doping is extremely 
inhomogeneous, as it contains a large variety of energy 
scales which depend exponentially on the spatial distri- 
bution of the impurities. In this paper we study the 
evolution of the OBD state upon application of a mag- 
netic field, which represents a straightforward experimen- 
tal probe of energy scales in magnetic systems. Precisely 
due to the large distribution of the effective couplings 
between the FMs, we find an amazingly rich response of 
the system to the applied field. The field scan reveals 
that the long-range order in the system is extremely ten- 
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FIG. 1: (Color online) Quantum phase transition between 
the order-by-disorder (OBD) phase and the novel disordered- 
free-moment (DFM) phase in a diluted coupled-dimer sys- 
tem. Left panel: At zero applied field, spins on intact dimers 
form singlets (solid ellipses), dimers of free moments (FMs) 
have a strong singlet component (dashed ellipse), whereas the 
other FMs (blue arrows) participate in the OBD state. Right 
panel: Upon applying a field, the OBD-FMs are mostly polar- 
ized, but local singlets and localized down-spins can survive 
on clustered FMs, leading to the DFM phase. 

uous for large spin gaps, and its field-driven destruction 
leaves behind local singlets (or spins oppositely polar- 
ized to the field) on even- (or odd-)numbered clusters of 
FMs, which are coupled at energies higher than the scale 
characteristic for Neel order (see Fig. [TJ. The local po- 
larization fields for these FM clusters cover a continuous 
range, so that the magnetization process with increasing 
field continues also after destruction of the OBD, and 
the resulting disordered-free-moment (DFM) phase (l3j 
is gapless. The DFM phase persists up to the sizable 
field which polarizes all the FMs, and the magnetization 
curve within this phase shows prominent features of in- 
termediate psewrfo-plateaus, still retaining a finite albeit 
extremely small slope, related to the distribution of the 
strongly interacting clusters of FMs. 

To quantitatively investigate the field response of a 
site-diluted spin-gapped antiferromagnet, we focus our 
attention on a two-dimensional 5 = 1/2 model of weakly 
coupled dimers 0, [H| , whose Hamiltonian reads 
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Here j runs over the two sublattices (A and £?) of a 
square lattice, x and y are the two lattice vectors, and Ci 
is the random dilution variable taking values and 1 with 
probability p and 1 — p respectively, h = g/isH is the 
applied field. The couplings J > J' determine the subset 
of strong antiferromagnetic bonds: for J'/ J < 0.523.. 
the bond anisotropy stabilizes a dimer-singlet ground 
state against the conventional Neel ordered state of the 
square-lattice antiferromagnet. All the results presented 
here refer to the field and doping effects deep within the 
dimer-singlet regime at J' / J = 1/4. 

The presence of lattice vacancies induces FMs which 
are localized in the vicinity of the unpaired spins 
which have lost their J-neighbor. Perturbation the- 
ory [13, O, E3] provides an effective coupling Jy « 
(— J\/r) cxp(— r/£o) between these FMs, where 
r — \i — j\, £o is the correlation length of the undoped 
system. We choose J\ w J'exp(l/£ ) in order for to 
correctly reproduce the limit J' for neighboring unpaired 
spins. For a deeper understanding of the Hamiltonian Eq. 
((TJ) , it is illuminating [l8| to study an effective model for 
the network of FMs, consisting of randomly distributed 
5 = 1/2 spins with effective couplings Jij 
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We investigate the original and the effective Hamil- 
tonian, given by Eq. (QJ and ((2]), using Stochastic Se- 
ries Expansion (SSE) Quantum Monte Carlo simulations 
based on the directed-loop algorithm \iSt\- For the orig- 
inal Hamiltonian Eq. ([I]) we study L x L lattices up to 
L=40 with dilution p = 1/8, whereas for the effective 
model Eq. ^ we randomly distribute spins on the same 
lattice sizes with a density p equal to that of the vacan- 
cies in the original model. Disorder averaging is typically 
performed over w 300 realizations. The ground-state 
properties are systematically obtained using a /3-doubling 
approach [2(J. Inverse temperatures up to (3 J = 2 15 
are necessary to observe the physical T — + behavior. 
In the following, we focus our attention on the uniform 
magnetization per spin m u — 1/N S Ylii^i )> wnere N s 
is the total number of spins in the system considered; 
on the uniform susceptibility Xu/J = dm u /dh; on the 
staggered magnetization m s = y S^(tt, ir)/L 2 , where 
S ± (n, tt) = l/(2i 2 ) EijiSfSf + SfS^} is the transverse 
static structure factor; on the correlation length £, ex- 
tracted from the g-dependent structure factor; and on 
the superfluid density p s = 1/ (2f3J)(W x + Wy), where 
W x / y -\ are the winding numbers of the SSE worldlincs. 
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FIG. 2: (Color online) Field scan of uniform magnetization, 
uniform susceptibility, and staggered magnetization in the di- 
luted system (a) and in the effective model (b). The lower 
right panel in (a) shows the temperature dependence of the 
uniform susceptibility at h — 0.175 J, indicating a saturation 
to a small but finite value as T — > 0. PL=plateau, BG=Bose 
glass, SF=superfluid. 

In the absence of a magnetic field, the effective cou- 
plings Jij give rise to long-range magnetic order of the 
network of coupled free moments. In particular, for 
J'/J=l/4 and p=l/8 we extrapolate a staggered mag- 
netization m s = 0.032(3) in the thermodynamic limit. 
Although the couplings Jij range between and J', the 
average coupling strength responsible for the long-range 
order turns out to be much smaller than J' [13]. Fig. 
[5] shows the evolution of the ordered moment under ap- 
plication of a field. It reveals that the antiferromagnetic 
order is already destroyed at a field h <C J' , namely at 
h = Iiufm = 0.007(1). Yet, a striking feature of this 
disordered phase is that the destruction of long-range or- 
der is not accompanied by the full polarization of the 
FMs, as it would ordinarily happen in a homogeneous 
antiferromagnet. At the critical field /idfm the uniform 
magnetization is found to be m u = 0.0208(6), much less 
than the value m u = pS = 1/16 corresponding to fully 
polarized FMs, which is attained at a much larger field 
^plateau/-/ ~ 0.5. Consequently the DFM phase, ap- 
pearing between /idfm and /ipiateauj is highly unconven- 
tional, retaining a finite uniform susceptibility and a gap- 
less spectrum. These unconventional properties have also 
been observed in the magnetic Bose-glass phase of triplet 



quasiparticles living on intact dimers [13|, [lj, UM, |21|. 
Later we will argue that these two phases bear indeed 
strong analogies, although involving different degrees of 
freedom. 

For h > /ipiatoau the saturation of the magnetization 
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of FMs leads to the full restoration of a gapped disor- 
dered phase due to the field 17]. Once the field reaches 
the value corresponding to the gap of the clean system, 
h = A = 0.60(1), a further Bose-glass phase is estab- 
lished in which rare clean regions develop a local mag- 
netization without the appearence of spontaneous order, 
corresponding to localized triplet quasiparticlcs 
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A derealization transition of the triplet bosons into a 
superfluid condensate corresponds to a further onset of 
lonR-range transverse order (m s > 0) at even higher fields 
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As shown in Fig. [2l the main features of the field 
dependence m u and m s for h < ft. p i a tcau in the doped 
coupled-dimer model of Eq. (fTJ) are very well reproduced 
by the effective model Eq. IJ3J, for which we take £o = 1 
as found by simulations at h — and p — 0. In particular, 
the fundamental appearence of a DFM phase with m s = 
and Xu > is confirmed in the effective model. This 
reveals that the FMs are essentially the only degrees of 
freedom responding to a field h < /i p i a tcau in the doped 
coupled-dimer system. 
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FIG. 3: (Color online) Finite-size scaling of the correlation 
length, staggered magnetization and spin stiffness in the vicin- 
ity of the quantum critical point hoFM- 

The novel quantum phase transition (QPT) be- 
tween the OBD and the DFM phase is studied us- 
ing finite-size scaling analysis of the correlation length, 
£ = LF^l^^Sh], of the superfluid density p s — 
L d ~ 2+z F Ps [L 1 /' J Sh], and of the staggered magnetization 
m s = L~P/ v F m3 [L^-^Sh], where Sh — h— /idfm, as shown 
in Fig. [3] This allows us to extract the critical expo- 
nents z, v and (3. For the original Hamiltonian Eq. ((T|) 
we find z = 2.0(1), v = 1.0(1) and = 0.9(1). These 
estimates are also confirmed in the effective FM model. 
The above exponents are fully consistent with those of 
the 2c? superfluid-to-Bose-glass (SF-BG) QPT previously 
studied in diluted bilayers and with the exponents 
found at the BG-SF transition for higher fields for the 
model Eq. (fTJ) 22]. In particular z is in agreement with 
the general theoretical prediction z = d [14| . and v sat- 
isfies the fundamental Harris criterion v > 2/d. Alto- 
gether the present results and those of Refs. [l3|, |22| point 
towards a general SF-BG universality in d = 2 for order- 



disorder transitions at which the uniform susceptibility 
Xu remains finite, corresponding to the absence of a gap. 

In the DFM phase, the magnetization curves of the 
original Hamiltonian Eq. (fl]) and the effective model Eq. 
@ show a dramatic feature: beside the large pleateau 
appearing at h > ft. p i a tcau, one observes the presence 
of apparent intermediate plateaus at around 3/4 and 
95% of the saturation magnetization. A detailed study 
of the temperature-dependent susceptibility in this field 
region reveals that these features are actually pseudo- 
plateaus (PPs), which retain an extremely small slope 
(Fig. HD. For both H and H FM the first PP extends 
up to h w 0.7 J'; a second PP markedly appears around 
h « 1.2 J' for H (it is rounded off for Hfm [22|); the true 
saturation plateau is only attained at h « 2 J'. These 
fundamental features can be understood within the pic- 
ture of strongly interacting clusters of FMs in the DFM 
phase. As shown in Fig. [T] the zero-field OBD phase is 
essentially inhomogeneous due to the random nature of 
the couplings. A majority of FMs are spaced from each 
other by an average distance (r) = p^ 1 ^, large in the 
small dilution limit, and interact via weak average cou- 
plings (Jcft) ~ pJ' [HI- However, fluctuations in the 
spatial distribution of the impurities also lead to small 
clusters of free moments located on neighboring sites, and 
thus interacting with much stronger couplings J'. If anti- 
ferromagnetically coupled in even-numbered clusters, the 
strongly interacting FMs participate only marginally in 
the OBD state of the system, and have a significant sin- 
glet component in their ground state wave function. This 
is directly revealed in a histogram of the bond energies 
(Fig. gfb)) E b = J b (S hb ■ S 3 ,b) where J b = J, J 1 and 
(1,6), (2,6) are the two neighboring lattice sites partic- 
ipating in the bond b. Beside the peak at E b ss — 3J/4, 
corresponding to singlets on intact dimers, a further peak 
at E b sa — 3 J'/4 is observed, corresponding to FM dimcr 
singlets, as well as a peak at E b — J'/2 corresponding 
to FM trimers. 

Applying a magnetic field h > ( J e ff) has clearly the ef- 
fect of destroying the long-range order of the FMs, but at 
the same time the FM singlets are left intact while odd- 
numbered clusters are not fully polarized (Fig. [1]), with 
the fundamental consequence that the antiferromagnetic 
order disappears but the FMs are far from saturation. 
This is clearly seen in the histogram of the local magnetic 
moments (S?) of the unpaired spins only (Fig. |U(b)): for 
small fields, a double-peak structure appears with a peak 
at (S?) — S corresponding to fully polarized FMs, and 
a strong quantum peak at (Sf) = corresponding to 
FM singlets. The large tails for < (Sf) < S and for 
(Sf) < come instead from FMs in odd- numbered clus- 
ters. In fact, we can resolve two more peaks at (Sf) = 1/3 
and (Sf) = —1/6 at larger fields, which correspond to 
partially polarized spins in FM trimers. Local FM clus- 
ters have widely different local gaps to full polarization, 
both due to their geometric structure (dimers, trimers, 
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FIG. 4: (Color online) Histogram of bond energies (a) and 
local magnetizations for unpaired spins (b). 

quadrumers, etc.), and to the local field they experience 
from the other FMs. Yet the distribution of rare FM 
clusters clearly assigns dominant statistical weight to the 
dimers, and this simple geometric fact is the reason for 
the appearence of the first PP: the magnetization pro- 
cess nearly stops until the local gap of the dominant FM 
dimers is overcome at h < J'. Nonetheless, the slope re- 
mains finite because the FM dimers have a distribution 
of local gaps. The magnetization value and the field loca- 
tion of the PP can be quantitatively related to the statis- 
tics of FMs clustered in dimers [22j. Analogously, one 
can quantitatively associate the second plateau with the 
statistics of the FM trimers [22|. Higher-order plateaus 
associated with larger local polarization fields should be 
expected, but they cannot be resolved within the given 
numerical accuracy. 

From the above data a clear picture of the DFM phase 
emerges: in this phase, a majority of the FMs are polar- 
ized, but antiparallel spins exist, corresponding to rare 
FM clusters. Upon a spin-to-hardcore-boson transforma- 
tion, these antiparalell spins take the nature of bosonic 
spin-down quasiparticles (|-QPs) localized on rare re- 
gions of the lattice [22|. This as pec t connects with the 
ordinary picture of a Bose glass [14J, and it further en- 
dows the OBD-to-DFM quantum phase transition with 
the nature of a localization transition: in the OBD phase 
the J,-QPs form a superfluid condensate, which is pro- 
gressively depleted by the applied field (acting as a nega- 
tive chemical potential), up to the point where the j-QPs 
undergo localization into a Bose-glass state, losing super- 



fluidity but retaining compressibility, which corresponds 
to a finite uniform susceptibility. 

It is evident from the above results that the DFM phase 
is relevant for experiments on site-diluted spin-gapped 
antifcrromagnets which display an OBD phase in zero 
field. The fundamental condition for the observation of 
the DFM phase is that the spin gap of the pure system 
be much larger than the maximum energy scale of the 
FM interaction (average inter-dimer coupling in weakly 
coupled dimer systems, inter-chain coupling in Haldane 
chains). This condition is necessary to ensure that the 
physics of the field response of the FMs is well separated 
in energy from that of the field-induced ordered state (for 
a detailed discussion, see Ref. [l3h . Joint magnetometry 
and neutron scattering measurements at relatively low 
fields should be sufficient to fully pinpoint this phase by 
demonstrating the absence of spontaneous order and the 
finite susceptibility down to zero temperature. Further- 
more, NMR measurements can show the rich structure of 
the distribution of local magnetic moments, similar to the 
histogram of Fig. [4] The low-field location of the DFM 
phase and its strong physical signatures in the magnetic 
observables make it the most accessible novel disordered 
phase in quantum magnets with lattice randomness. 
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